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I. INTRODUCTION 

The quest for strategies for combatting decoherence is 
of paramount importance to the control of open quantum 
systems, particularly for quantum information operations 
|l|. A prevailing unitary strategy aimed at suppress- 
ing decoherence is dynamical decoupling (DD) @,|1, [1], 
which consists, in the case of a qubit, in the application 
of strong and fast pulses alternating along orthogonal 
Bloch-sphere axes, e.g., X and Z. In the frequency do- 
main, where the decoherence rate can be described as 
overlap between the spectra of the pulse-driven (mod- 
ulated) system and the bath [Bj], DD is tantamount to 
shifting the driven-system resonances beyond the bath 
cutoff frequencies. The DD efficacy can be enhanced for 
certain bath spectra upon choosing the timings of the 
pulses so as to reduce the low-frequency parts in the sys- 
tem spectrum and thus its overlap with the low-frequency 
portion of the bath spectrum [J]. DD sequences are in- 
herently binary, i.e., their pulsed control parameters are 
discretely switched on or off. Realistically, the finiteness 
of pulse durations and spacings sets an upper limit on 
the speed and fidelity of DD-assisted quantum gate op- 
erations @,HI1]. 

An alternative strategy formulated here in full gener- 
ality is analog unitary control of multidimensional sys- 
tems subject to any noise or decoherence. It is effected 
by a system Hamiltonian whose time-dependence is vari- 
ationally tailored to optimally perform a desired gate 
operation. The vast additional freedom of non-discrete 
(smooth) Hamiltonian parametrization significantly en- 
hances the efficacy of decoherence control under realistic 
constraints compatible with the non-Markov time scales 
required for such control. Its formulation meets the long- 
standing conceptual challenge of simultaneously control- 
ling non-commuting system operators subject to noise 
along orthogonal axes. This is here achieved by working 
in an optimally rotated, different basis at each instant. 
The price we pay for such general optimal control is the 
need for at least partial knowledge of the bath or noise 
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spectrum, which is experimentally accessible @ without 
the need for microscopic models. The goal is to minimize 
its overlap with the spectrum of the controlled system, 
as was already shown for pure dephasing of qubits [3]. 



II. GATE ERROR 

We assume that the system Hamiltonian H$(t) imple- 
ments a desired quantum gate operation at time t, and 
aim at designing it so as to minimize the decoherence 
and noise errors. The system-bath interaction Hi then 
acquires time-dependence in the interaction picture un- 
der the action of H$(t) and the bath Hamiltonian H^. 
Assuming factorized initial states of the system and the 
bath, gtot(0) = g(0) <£> Qb, tracing over the bath, and fur- 
ther assuming that TrB[#i(i)£B] =0, yields for the sys- 
tem state g(t) the integrated (exact) deviation from the 
initial state (App.A), 

g(t) = (3(0) - AQ(t), 

Ag(t) = [dh [ di 2 Tr B [i?i(ti),[ffife),etot(i2)]]. (1) 
Jo Jo 

In what follows, we assume that up to t, the combined 
system-bath state changes only weakly compared to Hi, 
so that we approximate in (P) gtot(h) ~ £tot(0) in the 
integral. This means that the control is assumed effective 
enough to allow only small errors, consistently with the 
first order approximation of the solutions of both the 
Nakajima-Zwanzig and the time-convolutionless master 
equations H,Q. 

To justify this assumption, wc try to reduce the dis- 
crepancy between the states evolved for time t in the 
presence and absence of the bath by minimizing (Ag(t)) 
= (M/ 1 Ag(f)l^) averaged over all initial states \^) that 
are unknown in general. For a d- level system this av- 
eraging is tantamount to taking the expectation value 
with respect to the maximum entropy state g§ = dr 1 !. 
Assuming that Trs[#i(i)gs] =0, wc obtain our measure 
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of dccoherence (error) in the form of (App.A) 

(Ag(t)) = 2nRe f dh f ' dt 2 (Hi(ti)Hi(h)) 
Jo Jo 

rt - 1 2 \ 

J / SB 



Alternatively, we can rewrite ([8]) as 



SB 



(2) 



where n = l-(d+l) \ (-) SB = Tr SB [(-)£s ® Qb]- Hence, 
(Ag(t)) is always positive and proportional to the mean 
square of the interaction energy as observed in the inter- 
action picture (by a co-rotating observer). 

Since our aim is to suppress (Ag(t)} by system manip- 
ulations alone, we now separate system and bath parts 
by decomposing any interaction Hamiltonian in an or- 
thogonal basis of system states \j) as 



d 2 -i 

3=1 



(3) 



where the Hermitian Bj and Sj are bath and system op- 
erators, respectively, assumed to obey (Bj(t))B =TrSj(i) 
= and carry no explicit time dependence. In the inter- 
action picture 



= c iHBt BiC~ iHBt , 



Bj(t) 

U(t) = T+e-W 6 *W 



(4) 



We shall minimize (Ag(t)) for given, experimentally ac- 
cessible [(|, bath correlations 

*,•*(<) = (4(t)B fc > B - (5) 
It is expedient to define the decoherence matrix 



R(t 1 ,t 2 ) = e T (t 1 )Q(t 1 -t 2 )e(t 2 ), 



(6) 



which obeys R'(ti, t 2 ) = R(t 2 , t\). It is the matrix prod- 
uct of the bath correlation matrix J? formed from the 
coefficients $jk in ([5|) and the system-modulation (rota- 
tion) matrix defined as 



d z -x 
fc=i 

e jk (t) = (7) 



where we have assumed that r £r(SjSk) = 28jif The trans- 
formation ([7]) is at the heart of the treatment: it defines 
the instantaneous rotating frame where the system and 
bath are maximally decoupled, as shown below. 
We can now write ([2]) as (App.B) 



{Ag(t))=2^ / dti / dt 2 TrR(h, 
" Jo Jo 



(Afl(t))=44- 



dwTr[G(w)F t (w)], 



(9) 



t 2 ). 



(8) 



i.e., as the spectral overlap of two matrix- valued func- 
tions: the bath coupling spectral matrix G(u>) = 
dt c luJt RcJ?(i), and the system-modulation spectral 
matrix at finite time t [cf. ©] F t (ui) = je t (oj)et' (oj), 
e t {uj) = f Q dr c IWT e(r). In ([9]) wc have made use of 

the fact that t) = J? (t), so that it is sufficient to 
integrate over positive frequencies. 

Equation ((9|) constitutes a generalization of the "uni- 
versal formula" [f| to arbitrary multidimensional systems 
and baths. It provides a major insight: the system and 
bath spectra (all matrix components) must be anticorre- 
lated, i.e., Gjk(io) minima must coincide with (F t )jk(uj) 
maxima and vice versa to minimize ([9]), as illustrated 
below. It should be emphasized that for given w, both 
G(u>) and F t (u>) are positive matrices. Nevertheless, cer- 
tain components Gjk(ui), (Ft)jk(u>) may be negative if d 
> 2 (i.e., not for qubits). This may allow us to 'destruc- 
tively interfere' their contributions, i.e., engineer "dark 
states" or "decoherence- free" subspaces [ll|. These 
prospects of our general scheme will be explored else- 
where. 



III. DECOHERENCE MINIMIZATION 

Our goal is to find a system Hamiltonian Hg(ti), < 
t\ <t, implementing a given unitary gate U(t) at a fixed 
time t according to ((?]). This requires minimizing @ 
or ijSJ), (Ag(t)) -^min, i.e., minimizing the bath-induced 
state error in the interaction picture under Hs(t). We 
may similarly account for the effects of modulation or 
control noise, in addition to bath noise (App.C). 

The major difficulty in minimizing JS]) using (H])-® 
is that ((4|) involves time-ordered integration for arbi- 
trary bath and control axes. To circumvent this dif- 
ficulty, we make use of U(t\) instead of Hs, and as- 
sume a paramctrization U[fi(ti), t±] in terms of a set 
of real parameters fi(t±), which may be combined to 
a vector f(ti). The number of parameters may vary, 
since the paramctrization does not have to be complete. 
The boundary values /(0) and f(t) should be such that 

U(ti = 0)=I and U(ti=t) is the desired gate. 

If a bath coupling spectrum G(w) vanishes (has cutoff) 
at any high frequency, the overlap ([9]) can be presumed 
arbitrarily small under sufficiently rapid modulation of 
the Hamiltonian, such that all components of F t (uj) are 
shifted beyond this cutoff, thus achieving DD~0, y, 0] • 
Yet this may require a diverging system energy Further- 
more, fidelity generally drops with modulation energy, as 
discussed below. We therefore impose an energy con- 
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straint on the modulated system 

E s = [ dti(fl|(fi)) s = const, (10) 
Jo 

where (-)s = Tr[(-)d _1 i] [cf. ©]. An alternative con- 
straint 

E= ( dti |/(ti)| 2 = const. (11) 
Jo 

allows a simplified treatment. In general, E accounts for 
the fact that the time dependence of a parametrization 
cannot be arbitrarily fast and hence bounds the modu- 
lated Hs(ti), thus also limiting E$. 

The minimization of ([5]) subject to (fTTj) is an extremal 
problem in terms of /. Denoting by 8 the total variation 
with respect to /, the stationary condition can be formu- 
lated in terms of a Lagrange multiplier A as S_(Ag(t)) + 
XS_E = 0. Then, using the parametrization in .R [Eq. flBJ], 
Ve= {^-e(ii)}, yields the Euler-Lagrange equation 

'fih) = Xg(h), g{tx) = f dt 2 VRcTrR(t u t 2 ), (12) 

Jo 

where A is related to the constraint (fTTj) on E (App.D). 

We conclude the general treatment by recapitulating 
on the steps to find the optimal modulation of Hs{t\): 
1) After defining the 'cycle time' t and gate operation 
U(t), we declare a parametrization U[f{t\), t{\ which in- 
duces a parametrization e[/(ii),ii] that in turn yields 
R(ti,t 2 ) as a functional of / via (J3I)-(J7J), using our knowl- 
edge of ([5]). 2) We now solve (fl2|) for a given initial 
f i it (ti) satisfying the boundary conditions, e.g., such 

that #[/ init (*i)»*i] = [#(*)] ^, and calculate (Ag(t)). 3) 
The optimization is repeated for different values of A and 
Es in (|TT7|) is calculated for each solution. Among all so- 
lutions for which (Ag(t)} falls below a desired threshold 
value, we choose the one corresponding to the lowest Es- 
4) The chosen solution f(t±) is inserted into U[f(t{), t\] 
in yielding the instantaneous control parameters 

H s (h) = Y^"i(ti)Sj, ujih) = ^TriSjHsit!)}. (13) 

3 

IV. APPLICATION TO A QUBIT 

To apply the general procedure to a qubit for which Sj 
= &j (j = x, y, z) in ()13|) . we resort to the Euler rotation- 
angle parametrization, 

U(t) = C _ 5/3(*)*3 C -5/2( t )*2 c -3/l(t)*3^ 

In (fl~3|) . u>3(t) is now the level splitting, whereas w 1 ( 2 )(i) 
are Rabi flipping rates. We choose two examples of un- 
corrected (i.e., diagonal) baths, namely, an Ohmic bath 
with different cutoffs in A, Y, Z , and a Lorentzian noise 



spectrum superposed with a second Lorentzian such that 
a spectral 'hole' is obtained at different frequencies in 
X, Y, and Z. The corresponding bath coupling spec- 
tra are shown in Fig. [TJ along with our optimized mod- 
ulation spectra, which are contrasted with Uhrig's DD 
pulse-sequence spectra j3] (App.E,F). 




FIG. 1: Spectral overlaps between bath spectra Gi(u>) (solid 
red), modulation spectra i ? i OPT (w) (dashed green) for an op- 
timized 7r-gate at Es = 133.2 [(a),(b),(c)], and an optimized 
identity (0-) gate at Es — 181.1 [(d),(e),(f)], respectively, 
and modulation spectra fP UDD f or p U l se sequences (App.F) 
CUDD3 [(a),(b),(c)] and CUDD2 [(d),(e),(f)] (dotted blue), 
with i = 1, 2, 3 corresponding to X, Y , and Z- component, re- 
spectively. Graphs (a),(b),(c) represent an Ohmic bath spec- 
trum with softened cutoff, whereas graphs (d),(e),(f) repre- 
sent a Lorentzian spectrum with a dip. The optimal mod- 
ulation spectra i ? i OPT (tj) are always anticorrelated with the 
bath spectra Gi(u>). By contrast, Uhrig's pulse sequence spec- 
tra F? VD ~ D are only anticorrelated with Gi for Ohmic baths 
(a),(b) but not for the bath spectra (d),(e),(f). 

The minimized gate error is shown in Fig. H as a 
function of the energy constraint (fTU|) for both baths. 
Its comparison with the gate error obtained using var- 
ious DD pulse sequences reveals two differences. The 
first concerns the energy scale: in rectangular DD pulse 
sequences, each 7r-pulsc of duration T contributes an 
amount 7r 2 /(4T) to (JTUJ), which diverges for ideal pulses, 
T^O. By contrast, our approach assumes finite, much 
smaller E$ ■ The second difference concerns energy mono- 
tonicity: DD-sequences are designed a priori, regardless 
of the bath-spectrum, and hence only significantly reduce 
the gate error if Es has risen above some threshold which 
is needed to shift all system frequencies beyond the bath 
cutoffs ji}. In contrast, our approach starts to reduce the 
gate error as soon as Es > 0, since it optimizes the use of 
the available energy, by anti-correlating the modulation 
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FIG. 2: Qubit gate error (Ag(t)) in units of a reference 
{Ag(t)) rc{ which corresponds to a time-independent initial 
system Hamiltonian (with which the optimization started), 
as a function of the constraint Es. Solid red (dashed green): 
optimized identity (7r)-gate, Solid blue: periodic X-Z-"bang 
bang" (2,4,. . .,30 pulses). Separate points: concatenated 
Uhrig and related pulse sequences (App.F). (a) and (b) cor- 
respond to bath spectra shown on the left and right in Fig. [T] 
A different scale of -Es is used for pulse sequences whose Es is 
given in units of E^ /10, where E^ =tt 2 /(4T) = 2.467- 10 :! is 
the single 7r-pulse energy, assuming nearly-ideal square pulses. 



and bath spectra. 

We next consider the gate fidelity limitations as a func- 
tion of Es posed by leakage [l^] to levels outside the 
relevant subspace (here a qubit). In a 3- level A-system, 
any off-resonant control field acting on the qubit levels 
|1), 1 2), causes leakage to the unwanted level |3) [HI]. 
Such leakage and the ensuing incoherent decay |3) — > |1) 
incur gate errors that grow with Es (App.G). This be- 
haviour is illustrated in Fig. [3J which reveals that leak- 
age error is the more dramatic, the more energetic the 
7r-pulse sequences are. If 7r-pulses are experimentally im- 
plemented as (2m + Impulses, m>10 3 , it can therefore 
be expected that isolated manipulation on a subspace is 
difficult. This, together with the qubit-gate optimiza- 
tion, the general expressions for the gate error ([5]) and 
([{?[) and its minimization (|I2p are the main results of this 
work. 
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FIG. 3: Relative surplus error incurred by the inevitable leak- 
age to the additional level |3) (inset) caused by control f(t) as 
a function of Es, i.e., error with allowance for leakage com- 
pared to the error without (disregarding) leakage (App.G). 
Solid red: optimized identity-gate, dashed blue: periodic X- 
Z-"bang bang" (2,4,. . .,30 pulses) as shown in Fig. [21 a) and 
a bath spectrum as shown in Fig. [TJa),(b),(c). A (truncated) 
1/u) bath coupling spectrum describes an amplitude coupling 
between levels |2) and |3) (leakage). As in Fig. [2] a different 
scale of -Es is used for the pulse sequences: Es is here given 
in units of (3/I0)E^ 1) . 



V. CONCLUSIONS 

A) We have expressed an arbitrary gate error for finite- 
dimensional quantum systems as the spectral overlap be- 
tween the driven-system and the bath spectra. B) We 
have derived a non-Markovian Euler-Lagrange equation 
for the time dependence of control parameters whose so- 
lution maximizes the gate fidelity C) This solution leads 
to anticorrelation of the system and bath spectra. Hence, 
while DD-based methods rely on shifting the entire spec- 
trum of the system beyond that of the bath, our optimiza- 
tion takes advantage of gaps or dips of the bath spectra. 
D) The treatment of a qubit demonstrates that our ap- 
proach is significantly more economic in terms of energy 
investment than DD-bascd methods. Such energy saving 
may be crucial in terms of fidelity as excessive energies 
lead to leakage into additional levels [l3| , or increase the 
control noise [||. 
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APPENDIX A: DERIVATION OF THE 
DECOHERENCE (ERROR) EXPRESSION © 

The von Neumann equation for the total density oper- 
ator of the bath and system combined in the interaction 
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picture. 



^&rt(t) = -i[-ffi(*). &<*(*)]. (Al) 



can be written in integrated form 

etot(t) = gtot(o) - i / dh[Hi(ti), etot(ti)]. (A2) 

Jo 

Substituting (TA"2)) back into (|ATj) gives 

a 



ftot(*) = -i[#i(i),&ot(0)] 



dti[ffi(O,[i?i(ti),0t°t(*i)]], (A3) 



and after tracing over the bath. 
d 



g{t) = - / dtiTr B [ffi(t), [ffi(ti), ftot(ti)]]. (A4) 



9t ! 



Although we do not make use of the differential equa- 
tion for the system state g(t), it may be useful to men- 
tion that it can be obtained from (|A4[) by neglecting the 
bath correlations, i.e., setting gtot(^i) ~ £>b ® g{ti), which 
yields the second-order Nakajima-Zwanzig equation 
Replacing gtot(^i) ~ Qb <£> g(t) instead yields the second- 
order time-convolutionless equation. For the averaging, 
we make use of 



^mm^= ^^f^ (A5) 

[l4| to write the covariance of two operators A and B as 



Cov(A,B) = ($\AB\-$) - ($\A\V)($\B\$) 

= -<!)<% (A6) 

Expressing now the double commutator in ([1]) as 

[Hi{h), [fii(t 2 ),M = ([#i(*i),#i(*2)£tot] + ft.o.), and 
applying (|A6|) . we obtain ([2]). 



APPENDIX B: DERIVATION OF THE 
SPECTRAL OVERLAP ERROR © 



while U) reads 

d 2 -l „t rti 

Ag(t) = V dh dt 2 
j,k=i Jo Jo 



x {Qjkih-t^Sjih), S k (t 2 )g(0)} + h.a.} (B2) 
= J2 / dti / dt 2 

x {i? jft (ii,t 2 )[5 j ,^(0)] + /i.a.}. 

We can define a decoherence operator 
d 2 -i 

R(ti,t 2 ) = S j {t 1 )<f> jk (t 1 -t 2 )S k (t 2 ), (B3) 
j,k=i 

which obeys RJ(ti, t 2 ) = R(t 2 , t\)- Assuming finite d, @ 
then becomes 



(Ag(t)) = 2-Rc / dfci / dt 2 TrR(t!,t 2 ) 



o jo 
t ft 



o Jo 



d^ / dt 2 TU2(ti,t 2 ). (B4) 



Alternatively, by defining the spectral counterparts of the 
ingredients of (IB3I) : 



(B5) 



S s (w) = -= / dre^SjCr), (B6) 
y 1-K Jo 



/oo 
dt e iwt Re$^(t), 
-oo 



(B7) 



Equation (|B4[) can be written as the following spectral 
overlap 



(Ag(t)) = 2- / V Tt[St(w)G 3 - fc (w),§ fc (u;)] 

d J ° i,fc=i 

= ^-g] ^ G jk (u)F kj (u>). (B8) 



APPENDIX C: MODULATION ERRORS 



The differential equation (| A4|) for the system state can 
be written as [see comments following (|A4[) ] 

9 d 2 -i f 

— g(t) = - V / dti{$ ifc (t-ti)[5 i (t),5 , fc (ti)^(t)]+h.a.}, 
j,fc=i" /o 

(Bl) 



Since, in practice, a modulation can be realized only 
with finite accuracy, it is important to consider the ef- 
fect of modulation errors. To do so, we add to Hs{t) a 
small random Hamiltonian H^(t) which acts on the sys- 
tem variables and repeat the previous analysis without 
H^(t) in the interaction picture. In addition, we now per- 
form an ensemble average (also denoted with an overbar) 
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over different realizations of Hjq(t). Neglecting system- 
atic errors, Hn(t) = 0, we can in analogy to §5§ define a 
correlation matrix ffr N (ti,t2) with elements 



(CI) 
(C2) 



which gives rise to a noise contribution 

R N (t u t 2 ) = e T (ti)& N (t 1 ,t 2 )e(t 2 ), (C3) 

that must be added to (O with e(t) defined as before. 
Assuming ^ k (t 2 ,h) = *^.(ti,ta), we have R Ni {h,t 2 ) 

= R N (t 2 ,h), and © now holds for (Ag(t)): the double 
ovcrbar means that (Ag(t)) is averaged over both the ini- 
tial states and the ensemble. This analysis accounts for 
modulation errors if we use a modified correlation func- 
tion containing both system-noise and bath contributions 
and refer to the ensemble only. 



APPENDIX D: EULER-LAGRANGE 
VARIATIONAL ANALYSIS 



The minimization of (JSJ) subject to fjlOjl constitutes the 
original (unsimplified) extremal problem in terms of f. 
The stationary condition corresponding to (|10j) . 



S(Ag(t)) + XSE S = 0, 



(Dl) 



with variations fixed at the boundaries, 8f{ti)\t-L=o,t = 0, 
yields an Euler-Lagrange-equation 



RcTr 



U(ti)VU\h)-\ I dt 2 VR(t 1 ,t 2 ) 
o 



0. (D2) 



Here V; = d / dfi(t\) and the double dots denote a second 
derivative with regard to t\. In order to obtain (|D2[) . we 
have applied in (4) the relation 

H s (h) =i&(t 1 )U\t 1 ). (D3) 
The Lagrange multiplier in (TT2")) can be shown to obey 



yfb 2 + a(E-c)-b 



where 



dti 



*i 2 

dt 2 g(t 2 ) , 



o 



b = 



t fit 

dti / dt 2 /(0) ■ g{t 2 ) 



o Jo 
i 



(D4) 

(D5) 

(D6) 
(D7) 



Note that for /(0) = we have 6 = c = and (|D4j) reduces 
to A = y/E/a. 



APPENDIX E: BLOCH EQUATION ANALYSIS 

The state evolution of a qubit can be formulated 
in terms of the Bloch vector r with components rj = 
Ti[a jg(t)], j = 1,2,3, as the equation of a "top" forced 
by time-dependent torque 



t = L._ r + L + ■ (r - r ). 



(El) 



Here the matrix function 
ft 



L = 4Re f dti{fi T (t, h) - [TvR{t, (E2) 

has been decomposed into its (anti)symmctric parts L_± 
= {L±L T )/2, while 



Ho = -L+ b, 



bj = Aim J d*iTr[o: i; R(Mi)], 



(E3) 
(E4) 



is the quasi-steady state under the chosen time- 
dependent control. The term L + -(r— rg) accounts for the 
dynamically- modified relaxation of at non-Markov 
time-dependent rates that are the eigenvalues of L_At), 
reverting to the standard (Markov) rates 1 /Tj in the limit 
of slow control. The term L_-r = A.u> x r, reflects a bath- 
induced energy shift 

Auij = 2Re J^ dtiTr^fl^ti)], (E5) 

since it represents a unitary evolution observed in the 
instantaneous interaction picture. The elements of the 
SO([3|) generator matrices &_j can be calculated from 



2(<TAfc 



2i 



(E6) 



The optimized instantaneous control parameters u>j(t) 
are obtained upon minimizing the departure of r from 
its initial value and following the procedure in the main 
text leading to (fT3"|) . The results are illustrated in Fig. 0] 
(see also Fig. [T] main text). 



APPENDIX F: COMPARISON WITH UHRIG'S 
DD-SEQUENCE 

In Figs.[T]and[2]of the main text we compare our results 
with the following DD sequences: 

a) Concatenated DD (CDD) [3] defined by 

Pn+i = p n Xp n Zp n Xp n Z (Fl) 

with pq = f T denoting free evolution over time 
r, where pf DD = (fXfZ) 2 recovers periodic DD 
(PDD). 
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FIG. 4: Optimized time dependence of the control parame- 
ters of the system Hamiltonian, the solid red, dashed green, 
and dotted blue line show ^1,2,3, respectively: (i) parameters 
referring to the graphs (a),(b),(c) in Fig.l; (ii) parameters 
referring to graphs (d),(e),(f) in Fig.l. 



b) Uhrig-DD (UDD) [| defined (for n pulses in Z) 
by 

PT D = ft-T n Zf Tn _ Tn _,Z- ■ ■ Zf T2 _ T1 Zf Tl (F2) 
with 

n =isin 2 [7rj/(2(n + l))], 

where 

rf DD = fZf 

recovers the spin echo (SE) and 



P-2 



UDD = {fzf) 2 



(F3) 

(F4) 
(F5) 



the CPMG-sequence. 



c) Combined CDD and UDD in concatenated UDD 
(CUDD) Q defined by concatenating according to 



Pn+1 = PnXp n X 

an m-pulse UDD sequence 



po = pI db 



(F6) 
(F7) 



for m times. 



The named basic sequences can be iterated, i.e., repeat- 
edly applied. 



APPENDIX G: LEAKAGE FROM A SUBSPACE 

We can adapt our formalism to the situation where the 
d-dimcnsional state space (to which the relevant quan- 
tum information is to be confined) is a subspace of a 
TV-dimensional system state space [12]. To do so, the av- 
eraging of the initial states | VP) is performed on the sub- 
space, for which |*) = P|*}, where P = J2 n =i |<Pn)(<Pn| 
is the associated projector. Applying (|A5[) to {B2| and 
defining a matrix r — T} with elements 



Tr(SiPS k ) Tr(SiPSkP) + Tr(SiP)Tr(S k P) 



d(d + l) 



generalizes © to 



(Gl) 



(Ag(t))= I d*i / di 2 Tr[ ; R(ti,t2)r], (G2) 
Jo Jo 

which recovers (J5J for N — d, where = gxi^ife = 2§5ifc. 
Equivalcntly, 

{Ag(t)) = 2 f dwTr[e i t (^)G rc (^)e i (a;)Rer 
-^MG^HetHlmr] 
= t I dcTr[G tot MF t r M], (G3) 



which replaces (9). While G ro (w) = /^dt e iwt Re$(t) 
is identical to G(u>) in (9), here we also need G lm (ui) 
— /YJxjdt c 1LUt ImJ?(£), or the combined G tot (u>) = 

J^dt c iwt *(t), whereas in ^ F[_ (w) = j^w)Tet} (w) 
replaces F t (a;) in (9). 

(IG2I) encompasses both the internal decoherence ef- 
fects within the system-subspace associated with P and 
leakage effects related to a population (Q) = Tr[g(t)Q] of 
the orthogonal complement Q = I — P, averaged over all 
initial states on P, 



Tv[g(t)Q] = Tr[A^)P] 

= f dh [ dt 2 Tr i 2 )rrj 
Jo Jo 

Tl = Tl^ being a matrix with elements 



(G4) 
(G5) 



(T 



L it 



Tr(5,P5,.Q) 



(G6) 



If leakage is disregarded in the procedure minimizing 
(Ag(t)), it is likely that a stronger system modulation 
increases the population of Q, giving rise to a significant 
surplus error. This is illustrated in Fig. 3, where opti- 
mal and PDD-modulations originally designed within a 
two-level model [as shown in Fig. 2(a)] are reconsidered 
for a two-level subspace of a three-level system. This is 
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done by replacing the Pauli matrices &i with the corre- 
sponding Gcll-Mann matrices 7^ , multiplying U (t) with 
e -it/78 ^ separate the levels, and adding to Hi a leak- 
age term 76-E>l- The latter gives rise to an additional 
bath correlation function $l, assuming here that it can 
be described by a 1/w-bath coupling spectrum. The total 
system space is hence spanned by the energy states |1), 
|2), and |3), the projector onto the relevant subspace is 
P = J2 n= i \n)(n\, whereas Q — |3)(3|, and the states |\E') 
used for averaging are arbitrary superpositions of 1 1) and 



|2). The time-independent / is a parameter that con- 
trols the coupling to the "leakage bath". It reflects the 
fact that the energy of the leakage level |3) induces a free 
evolution, which is shifted to high frequencies for suffi- 
ciently large /, when |3) is strongly energy-detuned from 
the other two levels, thus providing a "natural" dynamic 
decoupling of our l/w-coupling spectrum, and hence the 
vanishing of the surplus error induced by leakage, justi- 
fying the two-level system approximation. 
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